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Interference and tunneling are two signature quantum effects that are often per- 
ceived as the yin and yang of quantum mechanics: particle simultaneously propagating 
along several distinct classical paths versus particle penetrating through a classically 
inaccessible region via a single least-action path. Here we demonstrate that the Dirac 
quasiparticles in graphene provide a dramatic departure from this paradigm. We 
show that Zener tunneling in gapped bilayer graphene (BLG), which governs trans- 
port through p-n heterojunctions, exhibits common-path interference that takes place 
under the tunnel barrier. Due to a symmetry peculiar to the BLG bandstructure, in- 
terfering tunneling paths form 'conjugate pairs', giving rise to high-contrast oscillations 
in transmission as a function of the gate-tunable bandgap and other control parameters 
of the junction. The common-path interference is solely due to forward-propagating 
waves; in contrast to Fabry-Perot-type interference in resonant tunneling structures it 
does not rely on multiple backscattering. The oscillations manifest themselves in the 
junction I-V characteristic as N-shaped branches with negative differential conductiv- 
ity, enabling new high-speed active-circuit devices with architectures which are not 
available in electronic semiconductor devices. 

PACS numbers: 



Quantum tunneling through two or more barriers that 
are placed closely together is characterized by transmi- 
sion which is sharply peaked about certain energies. Such 
'resonant tunneling' effect arises because particles can re- 
flect between the barriers and resonate at particular en- 
ergies, allowing enhanced transmission through the bar- 
riers. This resonance phenomenon is similar to that tak- 
ing place in optical Fabry-Perot resonators. Resonant 
tunneling is particularly desirable in applications since 
it can give rise to negative differential resistance-current 
that goes down as voltage goes up-an interesting behav- 
ior that can be harnessed to form new devices 0, Q ■ 

Here we propose an entirely different approach to re- 
alize oscillatory tunneling, which involves Zener tunnel- 
ing of Dirac particles through a p-n junction in gapped 
BLG [3, 4], a new material with a unique combination 
of electronic properties, such as the field effect and the 
possibility to open a bandgap by using external gates [H- 
7]. Interband (Zener) tunneling plays a crucial role in 
materials with several bands of carriers Q. Unlike the 
conventional tunneling through a potential barrier, which 
is controlled by the barrier properties, Zener tunneling is 
governed by an externally applied electric field that pro- 
duces mixing of states in different bands. Strong enough 
fields can induce interband transitions from the valence 
band of p-type material to the conduction band of n-type 
material, giving rise to tunneling currents. In conven- 
tional semiconductors, the tunneling rate is a monotonic 
function of the applied field F and the bandgap E g , given 
by an exponential exp(— nm^^Eg' 2 /2Fh) (here m is an 
effective mass) 0, [l(| ■ In a sharp departure from this be- 
havior, we find that transmission through a p-n junction 
in BLG oscillates as a function of the bandgap and ex- 
ternal field. The oscillations have 100% contrast, with 



transmission vanishing at particular nodal values of con- 
trol parameters (see FigJT]). 
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FIG. 1: Zener tunneling in BLG in the uniform-field model. 
Interference of two least-action tunneling paths results in os- 
cillations, n — 1,2.... Shown is transmission at normal in- 
cidence, p y = 0, as a function of bandgap size, in units 

Ao = ((i ? A) 2 /2m) 1 ^ 3 (semi-log scale). Numerical results 
(red symbols), obtained by integrating Eq.©, agree with 
the WKB result, Eqs.((T]),([4]) (blue curve) in the entire range 
of A, large and small. Inset shows schematic setup of p-n 
junction: the bandgap E g — 2A, the linear barrier potential 
V(x) = —Fx [see Eq.([2|l]. and a pair of interfering tunneling 
paths. 

The interference effects discussed below involve only 
forward-propagating waves and a single barrier, which 
makes them very diffeent from the Fabry-Perot reso- 
nances that arise from interference of waves undergoing 
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multiple reflection between several barriers. The origin 
of the oscillatory behavior can be elucidated by a semi- 
classical analysis of the dynamics in the barrier region. 
In contrast to the standard case of tunneling through a 
one-dimensional barrier, where a unique saddle-point tra- 
jectory in a classically forbidden region is found for each 
energy, here we obtain multiple trajectories. Further, the 
trajectories form pairs with complex conjugate WKB ac- 
tion values S and S*. Such pairs exhibit under-barrier 
interference, giving rise to an oscillatory transmission 
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where S — S' + iS". Both S' and S" are monotonic 
functions of the bandgap and field strength (see Eq.(0}). 
These oscillations will manifest themselves through neg- 
ative differential conductivity in the I-V characteristic. 

Relativistic-like behavior of carriers in graphene leads 
to many interesting transport phenomena 1 1 1!4 1 5j . How- 
ever, the oscillatory tunneling discussed here has not 
been anticipated by theory and is entirely different from 
Zener tunneling that governs transport in p-n junctions 
in semiconducting carbon nanotubes [l6l Il7| and sin- 
gle layer graphene [18j|. Theory of these systems (lil. l2pj 
yields exponential dependences that match closely those 
of Refs. 0, [T3| • Similar exponential dependence arises in 
the theory of p-n junctions in gated gapless graphene 
sheets 21|, with a momentum component along the p-n 
interface playing the role of a bandgap. 

The oscillatory tunneling in BLG opens door for de- 
signing new device arcitectures. Because the negative 
dl/dV arises solely due to single-particle effects, it is 
completely insensitive to the behavior in the doped re- 
gion. This represents a distinct advantage compared to 
resonant-tunneling (Esaki) diodes where the effect of 
temperature on dopants limit thermal stability. Also, the 
absence of multiple reflection makes the response time 
potentially much higher than for resonant tunneling in 
conventional quantum well heterostructures[l|. 

To clarify the origin of Eq. ([l} , we first consider trans- 
mission using the WKB formalism. Gapped BLG in the 
presence of a barrier potential V(x) is described by a 2 x 2 
quadratic Dirac Hamiltonian Q 
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where E g is the bandgap. We seek the wave function in 
the barrier region in the form ip(x) cx e h f<°o p ( x ^ dx ^, 
where ^ is a two-component spinor. The x dependence 
of momentum can be found from the energy integral E = 

± ( fp /2m) + A 2 ) +V(x). In the barrier region, 
—A < V(x) — E < A, this gives four complex roots 



Pj 



(x) = ±\J-P 2 y ± 2miy / A 2 - (V(x) - E) 2 , (3) 



where p y is a conserved y component of momentum. Two 
of the roots ([3]) have Imp > 0, while the other two have 
Imp < 0. Positive (negative) Imp correspond to decay- 
ing (growing) exponentials which describe particle prop- 
agation to the right and to the left, respectively. 

Focusing on the uniform-field model V(x) = —Fx (see 
FigfT] inset) and for simplicity setting p y = 0, we select 
from (|3|) the right-propagating solutions: p±(x) = (i ± 

l)?™ 1 / 2 (A 2 — (Fx — E) 2 ") 1 ^ 4 . These two solutions give 
complex conjugate WKB transition amplitudes e~ s / h , 
e -S l n where 



S, S* = {l±i)am 1 / 2 A 3 / 2 /F 



(4) 



with the prefactor expressed through the Euler beta func- 
tion, a = B(|,|) ~ !- 75 - 

The total transmission amplitude in the WKB approx- 
imation is the sum of the transmission amplitudes for the 
two decaying exponentials. Combining the contributions 
of the trajectories p±(x) we can write the WKB wave- 
function in the barrier region as a sum ae R *° + 
a*e B fxo p ~( x ^ dx . Interference between these evanescent 
solutions produces an oscillatory transmission amplitude 
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Since Re S — Re S* and Im S = — Im S*, the two contri- 
butions to the transmission amplitude are of equal mag- 
nitude and differ in phase by A8 = 2(j-lmS — ip). Here, 
ip = arg(a) is a phase offset between the two decay- 
ing exponentials which can in principle be obtained by 
matching solutions at the classical turning points, but 
in practice is more easily obtained through a numerical 
procedure, which gives <p s» 7r/2 (see below). 

For certain nodal values of the field strength F and 
the gap A the interference is destructive, and the trans- 
mission probability vanishes. Plugging the values (f?]) in 
Eq.©, we see that the transmission probability T = \A\ 2 
oscillates, going through nodes when am 1 / 2 A 3 / 2 /Ffi, is 
an integer multiple of 7r. This gives the nodal values 
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(Trn/a) 2 / 3 (F 2 h 2 /m) 1 / 3 , n = 1,2,3... (6) 



that match closely the nodes found numerically, FigfT] 

The oscillations in transmission, being a general fea- 
ture deriving from interference, are a robust and generic 
phenomenon. In particular, the oscillations do not re- 
quire a linear potential in the barrier region, and the 
WKB analysis may be straightforwardly generalized to 
an arbitrary potential profile V(x). Weak perturbations 
to the BLG dispersion also can be easily incorporated in 
the above analysis and shown not to matter as long as 
the perturbation strength is weak compared to the gap 
A. For example, the trigonal warping interaction can af- 
fect the dispersion within few meV of the Dirac point [3], 
thus its effect will be small in systems with gate-induced 
gap that can reach a few hundred meV 
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Another requirement on experimental systems in which 
the interference phenomena described above can be re- 
alized is that of ballistic transport in the p-n junction 
region. Recent observation of Fabry-Perot (FP) oscilla- 
tions in graphene p-n-p junctions [23j provides a clear 
signature of ballistic transport in this system. The os- 
cillation could be seen for the p-n interface separation of 
up to 60 nm, which sets a lower bound on the mean free 
path in the presence of a top gate. For a rough estimate, 
writing F = U/L with U a gate-induced potential dif- 
ference across a p-n junction and L the junction width 
(see Figj3] inset), from Eq.© we predict the number of 
experimentally accessible nodes 

am i/2 A 3/2 a A L h 

= -!==. (7) 
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For A = 100 meV, and using the effective mass in BLG 
m = 0.033 mo, we estimate the characteristic lengthscale 
£a ~ 3.18 nm. Taking eU = 4A and L = 60 nm, we 
arrive at n ~ 4, which indicates that oscillatory Zener 
tunneling is well within reach of current experiments. 

We now explain the origin of the oscillations from a 
different perspective, by mapping the transmission across 
the p-n junction to evolution of a two level system which 
is swept through an avoided level crossing. This alterna- 
tive formalism is specialized for the uniform-field model, 
and thus is less general than the WKB method. However, 
it provides intuition and affords an independent check on 
the WKB results by allowing us to numerically evaluate 
the transmission probability without any undetermined 
phase offsets. 

The key to this alternative formulation is an observa- 
tion that, for the uniform-field model V(x) = —Fx, the 
problem greatly simplifies in the momentum representa- 
tion. Indeed, since x = ifrd Px , the Schrodinger equation 
with the Hamiltonian turns into a first order differ- 
ential equation 
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where the cr, are the Pauli matrices in sublatticc 
space. This equation is identical to the time-dependent 
Schrodinger equation for a spin- 1/2 wavefunction with 
p x playing the role of time. 

There is a simple relation between the 'time evolution' 
governed by Eq.([8]) and interband transitions induced by 
Zener tunneling [24(. Asymptotically, at p x — ¥ ±oo, the 
eigenstates of the Hamiltonian are also eigenstates of a\ , 
having energies E ai= ±i = =fp x /2m. As we tune p x from 
-co to +oo, the system is swept through an avoided level 
crossing, as illustrated in Fig{2j Interband transitions are 
described by the process in which a state that started off 
in the o\ = — 1 cigenstate at p x = — oo will evolve into the 
<7i = +1 cigenstate at p x = +oo. The evolution is near- 
adiabatic at small F, with Zener tunneling described as 
(non-adiabatic) transitions across the gap. 
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FIG. 2: Evolution of a two-level system slowly driven through 
an avoided level crossing, Eq.([8j|. Non-adiabatic transitions 
between different levels, which correspond to Zener tunneling, 
take place in the Larmor precession region — pa < p < Pa, 
where pa = v / 2mA. Shown are adiabatic energy levels of the 
Hamiltonian, Eq.© (blue line) and schematic partition into 
regions of adiabatic evolution and Larmor precession. 



In this framework, the oscillations in transmission can 
be understood in a simple and intuitive way by noting 
that the Heisenberg evolution of momentum p x (t) cor- 
responds to sweeping through the avoided crossing at a 
constant speed, dp x /dt = F. Comparing different terms 
in Eq.([8|), we conclude that transitions may only hap- 
pen in the region — pa < Px J$ PA, where pa = v2mA 
(see Fig 12]), whereas outside this region the evolution is 
adiabatic (here we set p y — for simplicity). In the 
transition region the dominant term in the Hamiltonian 
is A<73. Spin rotation caused by this term can be de- 
scribed as Larmor precession about the z axis by an angle 
80 = (A/HF)pa- Periodic modulation of the transition 
rate of the form cos 56, resulting from Larmor precession, 
leads to an estimate of the oscillation period that agrees 
with the WKB result, Eqs.([T]),(|4]). 

The momentum- sweep analysis helps to understand 
the dramatic difference between transmission in bilayer 
junctions and single layer junctions. The latter problem 
can be mapped [25( to a canonical Landau Zener problem 
of a linear sweep through an avoided level crossing, for 
which transmission is a monotonic function of control pa- 
rameters exhibiting no oscillations. This is in agreement 
with the theory of p-n junctions in single-layer graphene 
21|. 

We now place this discussion on a firm quantitative 
ground by calculating the transition probability numer- 
ically. We solve the differential equation, Eq.®, in a 
suitably chosen interval p mm < Px < Pmax, taking as 
the initial state at p x = p m in the adiabatic ground state. 
From the numerical solution we determine the probability 
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to evolve into the excited state at p x — p max - The trans- 
mission probability, obtained in this manner for p y = 
and p m ax(min) = ±22pA, is shown in FigjTJ The results 
are compared with the prediction of the WKB approach, 
Eq.([T]), treating the prefactor \a\ 2 and the phase <p as fit- 
ting parameters. As illustrated in FigJTJ excellent agree- 
ment is found for the values ip = 1.6 and \a\ 2 = 0.78 
(which are tantalizingly close to ir/2 and 7r/4), indicat- 
ing that the WKB analysis provides reliable results. 

Integrating Eq.© at finite p y we find that the trans- 
mission oscillates and vanishes at nodal values of A in 
pretty much the same way as for zero p y . Comparing 
to the WKB analysis, which continues to apply at finite 
p y , we find that the WKB phase offset <p(jp y ) varies only 
weakly with p y . Using this numerical procedure, we may 
also straightforwardly take into account trigonal warping. 
Apart from a weak washing out of the nodes, we find no 
significant effect on the oscillations of transmission pro- 
vided the trigonal warping energy scale is less than the 
gap size. 
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FIG. 3: The I-V characteristic of a BLG p-n junction com- 
bines features of Esaki diode (N-shaped branches with nega- 
tive differential conductivity) and Zener diode (a breakdown- 
type behavior). Thus a single p-n junction can serve as an 
active circuit element with multiple functionality. Valleys in 
the I-V dependence correspond to n = 1 node of the oscilla- 
tions in transmission in Fig.l. Shown is the I-V dependence 
given by Eq. (|ll|l for parameter values: U/Vo = 0.1, 0.2, 0.3 
(curves 1, 2 and 3, respectively). Units are Uo = A(L/£&.) and 

la = 1Q' 4 N^{W/2tyIa)V , where W is the lateral width of 
the junction and N = 4 is the spin/valley degeneracy in BLG. 
Inset shows junction schematic, with U the built-in potential 
induced by doping or by gates, and E g — 2A the bandgap. 



Next, we proceed to show that the oscillatory tunneling 
reveals itself through distinct features in the I-V char- 
acteristic. The net tunneling current can be expressed, 
according to the Landauer formula [lij , as a sum of con- 
tributions of all conducting channels multiplied by energy 



distribution in reservoirs, giving 
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where W is the total length of the p-n interface, and the 
factor N = 4 is spin/ valley degeneracy in BLG. Here, 
accounting for the fact that transmission is dominated by 
small values of p y (see below), we treat the occupation 
numbers as p y independent and factor out the quantity 
T(F), the net transmission integrated over p y . 

Continuing to work with the uniform-field model, we 
treat transmission as energy independent and incorporate 
the source-drain voltage in the effective barrier potential 
via F = Fq + eV s d/ L (see Fig. 3 inset). Integrating over 
energies we have 



I = jV sd T(F + eV sd /L), F Q = eU/L. (11) 

The dependence of transmission T Py on p y may be found 
from Eq. (UJ) with S(p y ) and S*{p y ) evaluated using 
Eq.©. Since the transmission is exponentially small in 
the barrier width, and the width of the barrier region 
grows monotonically with p 2 , the net transmission T is 
dominated by small values of p y . Hence, we may approx- 
imate S and S* as 
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where a = a/2-B(|, |) « 2.4. Plugging these expressions 
in Eq.([T]) and performing Gaussian integration over p yi 
we find 
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<p = 2(p — f , where 5" and S" are given by Eq. (jH) . Based 
on numerical results, we ignored the p y dependence of the 
phase offset ip in Eq. (UJ) . Interestingly, the resulting I-V 
curve, Eq. (|ll|) . exhibits negative differential conductivity. 

A more accurate result for the net transmission T can 
be obtained by numerical integration of the exact WKB 
transmission over momenta (see Appendix). In that, the 
full dependence of S and S* on p y is retained, and also 
the contribution of the classically forbidden regions A < 



| Fx — E\ < yA 2 + (p 2 /2m) 2 is included, which is of 

subleading order in p 2 . 

The resulting I-V dependence is shown in Figj3] for 
several values of the 'built-in' (gate-induced) potential 
difference across p-n junction. Strikingly, the I-V char- 
acteristic combines features of the Zener diode (sharp 
rise of current above certain breakdown voltage) with re- 
shaped branches on which the differential conductivity is 
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negative, resembling the resonant-tunneling (Esaki) I-V 



characteristic 22J. Unlike the Esaki characteristic, the 



N-shaped branches occur simultaneously on the forward 
and reverse parts of the I-V dependence. The N-shaped 
features arise from oscillatory transmission (described 
by the uniform- field model), a mechanism very different 
from that leading to negative dl/dV in the Esaki diode. 
The valleys of current in Figl3] correspond to nodes of 
transmission (n — 1 in Fig.l). 

We note that p-n junctions of the type considered here 
can be realized using a configuration of gates which is 
already employed in current experiments [6, 23, 2(| 27 1. 
A minimal configuration is a dual-gate geometry with a 
wide back gate and a narrow top gate, such as that em- 
ployed in the work on FP oscillations 23] . Charging the 
two gates with voltages of opposite polarity, a bandgap 
can be induced under the top gate and, simultaneously, 
carrier density can be adjusted in the outer region. Ap- 
plying source-drain bias will produce l^d-dominated lat- 
eral electric field across the gapped region, corresponding 
to the regime U -C Vb where the effect of oscillations is 
most prominent (see Fig [3]). In addition, a built-in field 
field U can be induced by selective doping or by a third 
gate. 

In summary, transport in BLG p-n junctions is gov- 
erned by common-path interference under the tunnel bar- 
rier. Unlike Fabry-Perot interference that stems from 
multiple reflection between barriers, our interference ef- 
fect involves only forward-propagating paths and a sin- 
gle barrier. Common-path interference produces nodes in 
transmission as a function of the gate-tunable bandgap 
and other control parameters, leading to a complex I-V 
characteristic combining branches with negative differ- 
ential conductivity. The single-particle origin of nega- 
tive dl / dV makes it insensitive to the behavior in the 
doped regions which limits thermal stability and oper- 
ation speed of resonant tunneling (Esaki) diodes [22j . 
The operation speed is further enhanced compared to 
quantum-well-based devices by the absence of multiple 
backscatteringflj. We envision that BLG p-n junctions, 
owing to their multiple functionality and design simplic- 
ity, will become part of the future graphene electronics 
toolkit. 

We acknowledge useful discussions with N. Gu, C. M. 
Marcus and M. Rudner, and support from Naval Re- 
search Grant N00014-09-1-0724. 
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APPENDIX 

Our goal here is to calculate the net transmission T, 
a quantity used to evaluate the total current through p- 
n junction, Eqs.(9),(ll) of the main text. For that we 
evaluate transmission as a function of p y , where p y is the 
momentum parallel to the p-n junction. We will then 
integrate the transmission over p y to obtain the total 
integrated transmission through the p-n junction. Wc 
calculate this quantity by working in a WKB approxi- 
mation. 

When the WKB method is applied to the BLG p-n 
junction, the wave-function ip(x) is written as a sum 
of plane waves where the (potentially complex) wave- 
vectors are solutions of the classical equation 



h 2 ( K 2 {x) + k 2 y ) 

2m 



A 2 - (V(x) - E) 2 = (14) 



where k(x) is the WKB wavevector x component, E is 
total energy, hk y — p y is the momentum parallel to the 
p-n junction, 2A is the bandgap and V(x) is the gate 
potential. We are working in the uniform-field model 
V(x) = —Fx and, without loss of generality, set E = 0. 
The classical equation (|14[) has no real solutions for n 
in the classically forbidden region \x\ < Xf, where the 
turning points ±Xf(k y ) are determined by the condition 
n(±Xf) = 0. From Eq.(Hl]) we find 



A 

x f( k v) = -p\l l 



£a 



(15) 



where £ A is a lengthscale set by the gap. For \x\ < Xf, 
Eq. (fl4l has no real solutions and the wavefunction is en- 
tirely evanescent. 

In this problem, the forbidden region consists of two 
distinct parts. For \x\ < A/F, Eq. (fT4"|) has four com- 
plex solutions for k, given by Eq.(3) of the main text. 
Of these four complex solutions, two correspond to tun- 
neling from right to left and may be neglected, whereas 
the other two correspond to tunneling from left to right, 
with equal decay constants and a relative phase, which 
interfere when combined together. 

At | a: | = A/F, there is a doubly degenerate pure imag- 
inary solution to Eq. ([M)) , which corresponds to tunneling 
from left to right. In the outer part of the forbidden re- 
gion A/F < \x\ < Xf, the expression Eq. (fT4")l has four 



pure imaginary solutions, of which two correspond to 
tunneling from left to right. However, in this regime, the 
two tunneling paths are non-degenerate (have different 
decay constants, Re<Si ^ Re Si/), and we consider tun- 
neling only along the path with the longer decay length. 

In the semiclassical approximation, the amplitude of 
tunneling across the entire forbidden region, from x — 
—Xf to x — Xf is 



A(k y ) 



a exp — — 



-a exp — — 



exp 



2Si 

h 

(16) 

Here the action S is accumulated in the region —A/F< 
x < A/F. In this region there are two 'conjugate' tunnel- 
ing paths with actions S and S* , which interfere as dis- 
cussed in the main text. The action Si is accumulated in 
the outer regions —Xf < x < — A/F and A/F < x < Xf, 
where the tunneling paths are non-degenerate. In this 
region we consider only the tunneling path with longer 
decay length. The constants a and a* are parameters 
that may be found in principle by matching solutions at 
the classical turning points (but in practice are treated as 
fitting parameters with values obtained from numerical 
solution). The WKB actions are given by 



S 



Si 



A ' F dx 
-A/F ^a 
1 dx 

A/F ^A 



{Fx) 2 /A 2 , (17) 



p/2 
K y f -A 



^{Fx) 2 /A 2 -l (18) 



where £a and Xf are given by Eq.(fl5|). Note that the ac- 
tion S is complex (has real and imaginary parts) , whereas 
the action Si is pure real. The tunneling probabil- 
ity is given by the square of the tunneling amplitude, 
T{k y ) — \A{k y )\ 2 , and takes the form 



T(ky) 



f(k 2 y) = 



AAexp( -2— Ref(k 2 e 2 A 
A 



Ft A 

-Im f{ky(. Al 



F£ A ~ 

^k 2 y t\+i^fl-v 2 dv 



lfi(ky 



(19) 



(20) 



+2^ Fxf/A (k 2 y e 2 A - V^i) 1/2 dv 

where the first term describes the conribution of the re- 
gion \x\ < A/F, and the last term accounts for con- 
tributions of the outer regions —Xf < x < A/F and 
A/F < x < Xf (see Eq. (fT7j) ). Here Xf and £a are de- 
fined by Eq. (fT5|) . and the quantities <p{p y ) — 2 arg(a) and 
JV = 4|a| 2 are fitting parameters that may be found in 
principle by matching solutions at the classical turning 
points. Note that the tunneling probability is oscilla- 
tory at each k y , with the oscillations coming from the 
imaginary part of /. The imaginary part of / comes en- 
tirely from the region |a;| < A/F, where the action S is 
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complex, and where there are two degenerate tunneling 
paths, which interfere. 
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FIG. 4: Transmission (semi-log scale) as a function of p y and 
electric field. Units: pa = v2mA and Fo = A/£a- 

Fitting to the numerical solution of the BLG 
Hamiltonian, obtained by integrating time-dependent 
Schroedinger equation with a 2 x 2 Hamiltonian given 
by Eq.(8) of the main text, we obtain Af s» 3.2 and 
ip(k y = 0) « 1.6 (see Fig.l in the main text). 

The positions of the nodes, however, depend on k y . 
There are two main sources of this dependence. First, 
the period and phase of the oscillations is controled by 
the imaginary part Im/ dependence on k y . Second, the 
phase offset f(k y ) can vary with p y , producing additional 
shift of the nodes. However, since the contribution of 
Im/ to the net phase is greater than that of <p(k y ) by 
a large factor A/F£ A ^S> 1 (see Ed.([T9])). we expect the 
node positions variation with k y to be dominated by Im /. 

To compare the two effects, we use numerical solution 
to find transmission as a function of p y (see Fig.l). Map- 
ping out the nodes, we find that the phase offset tp varies 
only weakly with k y over the range ky£\ < 1 that dom- 
inates the integral. The position of the nodes is thus 
controlled mostly by modulation of the period. Hence, 
we approximate by taking ip sa 1.6 for all k y . 

We now wish to calculate 



. NW f 



(21) 



where we have assumed that the pn junction has lat- 
eral width W, and the factor of N = 4 arises from sum- 
ming over spins and valleys. We introduce the variable 
u = kyi\. The integral Eq. (l2~Tj) is dominated by small 
values of u. This is illustrated in Fig QJ which plots the 



transmission as a function of p y and electric field, as ob- 
tained from numerical calculation using the momentum 
sweep model. It is clear from the figure that transmission 
is dominated by small values of p y , i.e. u < 1. Therefore, 
we fit f(u) < 1 in Eq. (j2"0"f to a polynomial, and obtain 

f(u) fa 1.236-f-0.9u+0.45u 2 +i(1.236-0.8u+0.3u 2 )+O(u 3 ). 

(22) 

This second order polynomial fit provides an excellent 
approximation to f(u), as illustrated graphically in Fig [5] 




FIG. 5: Shows fit of function f(u), defined by Eq.pD). to the 
second order polynomial Eq. (|22[) . separately for the real and 
imaginary part of /. The solid red line is the exact curve, the 
dashed blue line is the polynomial fit. 

We then restrict the integral Eq.([2"Tj) to < U < 1, to 
obtain 



m NW 
T = 

2tt£ a 

where the function <!> is 

" du 



$(A/F£ A ) 



(23) 



<S>(A/F£ A ) = -^exp ( - 2— —(1.236 + 0.9m + 0.45z 
Jo V u V F£a 



' 2 f^(l. 



x cos" I 1.236- 0.8u + 0.3u 2 ) + 1.6 . (24) 



Evaluating numerically the integral over u and plug- 
ging the result in Eq. ([23"l) gives the total transmission, 
summed over p y . This result can now be used to obtain 
the I-V characteristic of the p-n junction, as discussed 
in the main text. 



